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We employ a quantum Langevin equation approach to establish Bloch-type dynamical equations, 
on a fully microscopic basis, to investigate the measurement of the state of a coupled quantum dot 
qubit by a nearby quantum point contact. The ensuing Bloch equations allow us to analyze qubit 
relaxation and decoherence induced by measurement, and also determine the noise spectrum of the 
meter output current with the help of a quantum regression theorem, at arbitrary bias-voltage and 
temperature. Our analyses provide a clear resolution of a recent debate concerning the occurrence 
of a quantum oscillation peak in the noise spectrum. 
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The quantum point contact (QPC) has been proposed 
and designed as an efficientjPractical weak continuous 
measurement device (meter) Recently, 
theoretical analyses of the detector noise spectrum have 
been carried out under arbitrary bias-voltage and tem- 
perature conditions using a perturbative Green's function 
methods, '0, also a quantum jump techniaue[lo|. and 
a Gurvitz-type master equation approach with a multi- 
subspace ansatz[ll|. in which some results have led to 
disagreements and remain inconsistent 1 1 2l| . Therefore, 
this problem requires further study. 

Here, we revisit this issue employing a recently devel- 
oped generic quantum Langevin equation approachj^^, 
which we approximate to establish a set of Bloch-type dy- 
namical equations describing the time evolutions of qubit 
variables explicitly in terms of the response and corre- 
lation functions of the detector variables. These Bloch 
equations provide analytical expressions for the bias- 
voltage and temperature T-dependent relaxation and de- 
coherence of the qubit induced by measurement. We then 
evaluate the detector current and frequency-independent 
shot noise within 2nd-order perturbation theory, and 
they are shown to fulfill the nonequilibrium fluctuation- 
dissipation (NFD) theorem[l3|. Finally, we calculate the 
frequency-dependent noise using a quantum regression 
theorem (QRT) based on the derived Bloch equations fisf. 

The system we examine is a double quantum dot 
charged qubit ci and C2 coupled with a nearby low- 
transparency QPC detector. In order to properly ac- 
count for dissipative effects in the qubit evolution, we 
write the total Hamiltonian of the system in terms of the 
eigenstate basis of the qubit, i.e. Ca — cos f ci 4- sin |-C2, 
c/3 = sin ^ci— cos -102, with d = tan~^ j {ft is the tunnel- 
coupling between two quantum dots with energies ±^): 
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where A = 2^/6^'+^ measures the energy difference 
between two eigen-energies, i.e. the Rabi frequency; cj^^. 
(cjjk) is the electron creation (annihilation) operator in 
lead ri{= L, R) of the QPG with a flat density of states 
P^; Ql = -Q^{t)cos^, Qt = = -iQ,(i)sini9, q\ = 
f Q,(i) with Qi(t) = xEk^qcLcflq and Q^it) = Q\{t) 
are the generalized co-ordinates; W and x are the di- 
rect and the qubit-modified tunneling matrix elements 
of the QPC, respectively, which are taken as constants; 
Fq^ = CTz = i(4c„ - clcp), Fqt = a± = 4(^)C^(„), 
and Fqi — 1 (the unit operator) are the generalized 
forces. Here, the terms Q^Fq± describe energy- exchange 
(inelastic) processes during measurement, in which a 
qubit-state transition takes place from one eigenstate 
to the other one due to the measurement. Note that 
our setup, Eo. ITll , is similar with those of previous 
studies[i[l|i[i3r 



In our derivation of Bloch-type dynamical equations 
for the qubit variables, operators of the qubit and the de- 
tector are first expressed formally by integration of their 
Heisenberg equations of motion (EOM) , exactly to all or- 
ders in the tunnel coupling and qubit-detector coupling 
constants, W and x- Next, under the assumption that 
the time scale of decay processes is much slower than that 
of free evolutions, we replace the time-dependent oper- 
ators involved in the integrals of these EOM's approxi- 
mately in terms of their free evolutions (Markov approx- 
imation). Finally, these EOM's are expanded in powers 
of coupling constants up to second order, taking the con- 
venient and compact form|l3l|: 
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where the summation notation c is over {z,±,l}, and 
n, / are over {1,2}; the statistical average {■ ■ ■)m{q) is 
in regard to the detector (qubit) variables; and the su- 
per(sub)script "o" signifies free evolution of the corre- 
sponding variables; t = t — t'. In conjunction with free 
evolution, the qubit dynamics are modified by the mea- 
surement process in a way that relates to the response 
function, Rni{t,t'), and correlation function, Cni{t,t'), 
of free meter variables, which are defined as: 



Ral{t,t') = ^e{T){[Qnoit),Qloit% 
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Cnl{t,t') -0{T){[Qno{t),Qlo{t')] + ),n- (4) 

The nonvanishing correlation functions can be readily ex- 
pressed in terms of reservoir Fermion distribution func- 
tions of the meter, and their Fourier transforms are|l.ll|: 
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-Rl2/2l(^) = gi{^^ ± V), Cx2l2\W) = 9\T'^ ^ 

(5) 

with gi = -KpLpiix^ /2 a-nd <p(x) = a; coth(a;/2). V is the 
bias-voltage applied between the left and right leads of 
the meter. We use units with h = ks = e = 1. 

For long time scale of interest, making the replacement 
/^^ dr =^ dr, these EOM's, Eq. ©, can be further 
simplified as [here, we write {(yz(±)) ^z(±)\'- 
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&± = ±iAa± - —a± +rd(7^ +T\\az + c\\, (6b) 
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laxation rate ^j- and decoherence rate ^ are given by 



2^, Fy = isin2?9C+(A), = -i sin^ i?i?+(A), c,, = 
isin2i9i?+(A) [C(i?)± = C(i?)i2 ± C(i?)2i], and the re- 



;i = sin2^C+(A), -l = -i-+cos2i?C+(0). (7) 
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It is readily seen that the relaxation time (Ti) stems 
completely from inelastic measurement events, which is 
conceptually consistent with the physical definition of 
qubit relaxation. These inelastic processes also con- 
tribute to qubit decoherence with the partial rate, l/2Ti. 
In contrast, elastic processes do not induce relaxation but 
do contribute to pure decoherence with the partial rate 
cos^ 'i?C+(0). In the case oi d — (no inter-dot hopping, 
= 0), the relaxation rate is, of course, zero, mean- 
ing that the qubit is completely localized, ^ — C+(0), 
and — (depending on the initial state), — 0; 
On the other hand, if -i? = 7r/2 [the symmetric (S) case, 
S — 0], ^ — C^{n) and ^ ~ We have the general 
steady-state solutions of Eqs. and (|6b|l as: 
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The Markov approximation employed in the derivation 
of Eqs. © requires that (1) measurement-induced decay 
of the qubit Td < A; and (2) rapid internal decay of the 
two electrodes, I/tc > V (coarse-graining assumption), 
which implies that the noise spectrum in the following 
calculation is meaningful only for low frequencies uj <V. 

The tunneling current operator through the meter is 
defined as the time rate of change of charge density, Nl — 
^j,c^j,CLk, in the left lead: 

JLit)^NL = -i[(Q? - Q^)<7, + (Q+ - Q+)(7_ 

+ {Q- -Q-)a+ + {Ql~Q'^)-i]. (9) 

Using linear-response theory we have 
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with Ai = \ cos^ i9 - ^ cos -daz 
In the case of = 0, the steady-state current is //( 
y{9o + ffi ± 2^goffi) with go = where 
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means the QPC current when dot 2(1) of the qubit is 
occupied by electron; For the S case S = 0, we have 



I^goV + giV 1 
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To address the noise spectrum of a "symmetric" detec- 
tor, we employ its definition as the Fourier transform of 
the current-current correlation function, S'(t), which can 
be calculated using linear-response theory. 



dre^-i([<5JL(t),Wi(t')]+)™,„ (12) 
1 ^ 



with 5JL{t) = Jbit) — {J Lit))- It is well-known that 
the shot noise spectrum consists of the w-independent 
Schottky noise (the pedestal) and the w-dependent 
part. Substituting the current operator of Eq. 10) into 
Eq. fl^ and using the definitions Eqs. ©-(0), we obtain 
the a;-indepcndent noise, S^)^ as: 
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where the first terms on the right hand side of Eq. H1U|) 
at t ^ oo and Eq. H13|l obey the NFD relation, while the 
other two terms represent the generalized NFD relation 
due to energy-exchange processes involved in the course 
of measurement Ti]. 

From Eq. (|10|l . it is clear that the time evolution 
of the qubit variables determines the temporal behav- 
ior of the meter current, and is consequently respon- 
sible for the frequency-relevant part of the noise spec- 
trum. To calculate the w-dependent noise, it is thus 
necessary to evaluate the two-time correlation functions 
(Tahir) = {(Tait + T)ab{t)) {a,b = z,p,m) using the QRT 
starting with the dynamic equations ®- The QRT states 
that given closed-form equations of motion of one-time 
averages of a set of operators Oj: 



J2LjAOkit + T))+\,, (14) 



then the two-time averages of Oj with any other operator, 
"P, also obey the same equations 15j: 
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{O, it + T)V{t)) - ^iMOk it + T)V{t)) + A, {V{t)). 

k 

(15) 

Because the solutions of interest here involve long time 
behavior after the initial turn-on, we replace the initial 
values of crafc(O) = {<Ja(t)ah(t)) and the last term {V{t)) in 
Eq. (|15|l with their respective steady-state values. Once 
these equations for aab{T) are solved and substituted 
into the definition, Ea. (|12|l . we calculate the frequency- 
dependent noise spectrum S(lu). 



There is no general analytic expression for the noise 
spectrum, except in two special cases: In the absence of 
coupling, n^O, lit) = [i+2(f )2]i?-(0)-f£a,(i), but 
(Tj = 0, therefore the total noise is independent of with 
S = Iii^r) coth(V^/2T); On the contrary, for the S system 
S = 0,we have /(<) = /q + /i(t) = (50 + gi)V - Ac7p{t) + 
Ba^it), with A ^ ^giV and B = C-{A). Moreover, 
since CT^(t) is decoupled from cr±(t) in Eqs. JHl due to 
Tz — and Fn — 0, we have o'pz(r) = (Tzp{t) = 0, and 
then S{t) = Re[A'^app{T) + B'^a^.^ij)] - {hf. Finally, 
the total noise spectrum, S{lu) ~ Sq + Si{uj) + S2{(^), is 
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with SI = AVy^gogi. The pedestal shot noise, 5*0, is 
identical with the previously found result 0. For the 
frequency- relevant noise spectrum in a S qubit, the de- 
coupled qubit dynamics give rise to two distinct compo- 
nents: Si{uj) stems from the transverse qubit dynamics, 
CTp, i.e. elastic and inelastic measurement-induced qubit 
decoherence, which is similar to the calculation of Ref. 
but augmented with our specific predictions for V- and 
T-dependent parameters, and SI; while S'2(w) is gen- 
erated by the longitudinal qubit dynamics Cz (qubit re- 
laxation), which is exclusively due to inelastic processes 
in measurement. Clearly, our Si{uj) is similar to Shnir- 
man's Eq. (35)0 except for the reduction factor. More- 
over, our S2{uj) coincides with C3(w) [Eq. (39)] under 
the same limit conditions, and also vanishes for F < A 
at T = 0. In our model, 5*2 (w) is at least one order of 
magnitude smaller than 5*1(0;), irrespective of T, and 
.91/30 [inset in Fig. 1(a)], indicating that the detector 
spectrum mainly refiects S qubit decoherence behavior. 
Noticing that our different low-bias behavior of 5*1(0;) as 
compared to Ref. is due to the coarse-graining assump- 
tion, which requires V ^ A to guarantee accurate noise 
spectrum at 0; ~ A. 

Figure 1 (a) plots measurement-induced qubit decoher- 
ence effects on the noise spectrum in the maximal meter- 
response case gi ~ go = 0.25. It may be observed that 
with increasing F^ (corresponding voltages V/A = 1.5, 
2.0, 3.0, and 4.0), S{uj) changes from a Lorentzian-type 
function, with a peak located at the qubit Rabi frequency 
o; = A, to the spectrum shape centered at w = 0, indicat- 
ing the Zeno effect regime. This behavior can be under- 
stood as follows: (1) a low voltage apphed to QPC can 
not provide sufficient dissipation (small F^) to suppress 
the Rabi oscillation of a S qubit, thus a peak occurs in 
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FIG. 1; (a) S{u!) as function of frequency u>/ A for various Td 
values. The parameters we use in calculation are: go = gi = 
0.25 and T = 0.1 A; (b) The peak-to-pedestal ratio vs. bias- 
voltage; (c) Ranges of various functional forms for the noise 
spectrum for a given QPC transparency go = 0.25. Inset: 
S{lu) vs. cj/A. Solid lines denote results for various gi/go 
values at go — 0.25. The dashed and dotted curves are for 
(?i — go but lower values of go = 0.125 and 0.06, respectively. 



the noise spectrum contrary to the previous prediction 
of Ref. as mentioned above; on the other hand, (2) a 
high voltage can freeze the qubit due to enhanced re- 
laxation, indicating appearance of the Zeno effect. The 
specific value of voltage dividing the two regimes nat- 
urally depends on the QPC-qubit coupling for a given 
QPC transparency go. Based on Eq. (|16b|) . the transi- 
tion occurs under the condition > A/V2, which leads 
to a rough transition boundary involving x/VF and V 
(log|^ oc —\ogV) as shown in Fig. 1(c) for the QPC 
with a fixed transparency ga — 0.25. The Lorentzian 
spectrum regime extends to a wider range of V for sys- 
tems with weaker detector-qubit coupling x- The inset 
of Fig. 1(c) exhibits the fact that the spectrum changes 
from the Zeno regime back to the Lorentzian regime as gi 
decreases. We should notice that the transition bound- 
ary depends on go due to gi < go, though is inde- 
pendent of gQ. In the high- voltage limit, V ^ A, T, we 
can still obtain the Lorenzian noise spectrum of detector 
output via applying a QPC with very low conductance, 
go <C A/y, to guarantee F^ ^ A. This remark is shown 
by dashed and dotted lines in the inset of Fig. 1(c). 

Another interesting feature is the quantum upper 
bound of 4 for signal-to-noise ratio (SNR) in the qubit 
measurement, which was first pointed out by Korotkov 
and Averin in the high voltage limit From Fig. 1(b), 
we observe that (1) the SNR increases rapidly initially 
with rising voltage, and reaches a maximum value less 
than 4 (depending on gi) around ^ = A; (2) increasing 
qubit-QPC coupling 511 results in reduction of the SNR 
due to enhancement of 5*0, and the SNR becomes 2 in 
the high voltage limit in the case of maximal response 
91 = go- These results indicate that low QPC trans- 



parency go, weak meter-qubit coupling, gi, and > A 
are necessary for an efficient meter. 

The efficiency of a meter also depends decisively on 
temperature*^. The spectrum peak at w = A in Fig. 1(a) 
will gradually disappear with increasing T (not shown 
here) , which is ascribed to enhanced detector-induced de- 
cay, Td- From Eq. 10, F^ is independent of F at T < A/2 
and V < A, implying that the QPC is always a good me- 
ter for a S Qubit at low temperature and bias-voltage. 




FIG. 2: (a) S{ij) vs. uj/A for an asymmetric qubit with 
S — Q at various bias-voltages; (b) S{uj) for different qubits 
at 1/ = l.OA. Other parameters: go = gi = 0.25, T = O.IA. 

The noise spectra for asymmetric (A) qubits are sum- 
marized in Fig. 2. In contrast to the S qubit, the A qubit 
longitudinal dynamics start to play a dominant role in 
the spectrum around uj — Q when V > A [Szz in the in- 
set of Fig. 2(a)], making S{lo) approach a zero-frequency 
maximum more quickly, albeit that the transverse spec- 
trum, Spp, still has a peak at w = A. This is to say that 
breakdown of the resonance condition causes the A qubit 
measurement to enter into the Zeno effect regime more 
quickly, even with F^ <C A. 

In conclusion, we have analyzed the noise output spec- 
trum of a QPC-qubit measurement system for arbi- 
trary voltage and temperature conditions using a gen- 
eral quantum-operator Langevin equation approach to 
derive Bloch equations for the system. In contrast to 
the results in Ref. Hol which improperly ignores quan- 
tum interference between elastic and inelastic processes 
of measurement, our calculations indicate that qubit os- 
cillations always cause a peak in the QPC noise spectrum 
at the Rabi frequency. A, for a S qubit, because the qubit 
coherence dynamics dominate contributions to the noise 
spectrum (if the condition F^ < A/y/2 is satisfied), which 
agrees with previous studies in the high voltage limit 0. 
However, the coherent peak at the Rabi frequency is sup- 
pressed by the relaxation dynamics of an A qubit even 
when Yd <SC A. Our analyses provide physical insight into 
qualitative criteria for design of an efficient detector: (1) 
low transparency go of a QPC; (2) weak meter-qubit cou- 
pling, gi, (3) small asymmetry ratio, 5/^', (4) relatively 
low temperature. 
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